Intrinsic spin Hall effect in platinum metal 
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Spin Hall effect (SHE) is studied with first-principles relativistic band calculations for platinum, 
which is one of the most important materials for metallic SHE and spintronics. We find that intrinsic 
spin Hall conductivity (SHC) is as large as ~ 2QQQ[h/ e)[Q,cm)~^ at low temperature, and decreases 
down to ^ 200(/i/e)(r2cm)~^ at room temperature. It is due to the resonant contribution from the 
spin-orbit splitting of the doubly degenerated d-bands at high-symmetry L and X points near the 
Fermi level. By modeling these near-degeneracies by an effective Hamiltonian, we show that SHC 
has a peak near the Fermi energy and that the vertex correction due to impurity scattering vanishes. 
We therefore argue that the large SHE observed experimentally in platinum is of intrinsic nature. 

PACS numbers: 71.15.Rf, 72.15.Eb, 72.25.Ba, 75.47.-m 



Spin Hall effect (SHE), i.e., the transverse spin cur- 
rent generation by the electric field, is an issue of in- 
tensive current interests both theoretically and experi- 
mentally since the theoretical proposal for its intrinsic 
mechanism in semiconductors [l|, I4I . This effect enables 
us to control spins without magnetic field or magnetic 
materials, which is a crucial step for spintronics. In ad- 
dition to semiconductors, the SHE in metallic systems is 
currently attracting interest, stimulated by experiments 
on the SHE or inverse spin Hall effect (ISHE), i.e. the 
transverse voltage drop due to the spin current jl, 0, 0] . 
SHE/ISHE in metals has the following importance and 
advantages compared with that in semiconductors; (i) A 
contact with a ferromagnetic metal does not suffer from 
conductance mismatch [6| , and one can make use of the 
spin-polarized current supplied from it. Thus, techniques 
developed in metallic spintronics can be utilized, (ii) The 
spin Hall conductivity (SHC) is much larger than that in 
semiconductors. The value of SHC obtained in Ref. 
is 4 orders of magnitude larger than that in GaAs 
Naively this appears to be attributed to the large num- 
ber of carriers, whereas the band structure is important 
as we discuss below, (iii) The Fermi degeneracy temper- 
ature is much higher than room temperature, and hence 
quantum coherence is more robust against thermal agi- 
tations than in semiconductors. We note that the spin 
diffusion length is relatively small in metals, e.g. 10 nm 
in platinum (Pt) causing fast decay of the SHE sig- 
nal. However, it is not a crucial obstacle for observation 
and application, by designing the device as demonstrated 
in Ref. 0]. 

Compared with the recent experimental advances in 
metallic SHE, its theoretical understanding is still lack- 
ing and is urgent. Among metallic systems, Pt shows re- 
markably large SHE surviving even up to room tempera- 



ture [1, 01 , whereas aluminum and copper show relatively 
tiny SHE The SHC in Pt at room temperature is 240 
[h/ e){Q.cm)~\, ten times larger than that of aluminum 
at 4.2K. In [J| this difference is attributed to a magni- 
tude of spin-orbit coupling for each metal. However, Pt 
seems to be special even among heavy elements, and the 
SHC does not simply scale with the size of the spin-orbit 
coupling. Such behavior cannot be explained within the 
extrinsic mechanism [loj , where material properties 
are represented by a few parameters such as the size of the 
spin-orbit coupling. This material dependence strongly 
suggests a crucial role of intrinsic contributions, which 
has been largely overlooked. It is thus highly desired to 
study the intrinsic SHE of Pt as a representative material 
for metallic SHE. This analysis opens up the possibility 
to theoretically design the SHE in metallic systems. 

This discussion on separating intrinsic and extrinsic 
mechanisms is analogous to the long-standing debates on 
the anomalous Hall effect (AHE) [ui, 0, 111 Q] ■ In semi- 
conductors, there have been exp erimental reports on the 
SHE in n-GaAs [Z], p-GaAs [I^ and n-type InGaN/GaN 
superlattices [l6| . It is now recognized that the SHE in n- 
type GaAs is due to the extrinsic mechanisms, i.e., skew 
scattering and side-jump contributio ns [9l.ll7l . while that 
in p-type GaAs is mostly intrinsic [iSl . Il9|. In metals, 
the conventional understanding has been that the skew 
scattering is dominant in AHE. However, recent studies 
have revealed that the intrinsic contribution can be dom- 
inant for AHE in metals when the Uxy is of the order of 
lO'^ (0cm) ~^ and the conductivity axx is in the range of 



lO^'-lO^O^^ cm- 



This dominant contribution of 



intrinsic mechanism is confirmed by the detailed compar- 
isons between the first-principles calculations 21, 22, 2§| 



and experiments 24 1. 



In this Letter we present an ab initio calculation for the 



2 



SHC in Pt, and its analysis based on an effective Hamil- 
tonian. We find that there are near-degeneracies near 
the Fermi level (Ep) at high-symmetry X and L points 
in the Brillouin zone (BZ) for the fee lattice. They give a 
prominent enhancement of SHC in Pt. We determine an 
effective Hamiltonian near X and L points, and demon- 
strate robustness of the SHE against impmities. 

The band structure of Pt is calculated using a fully rel- 
ativistic extensi on [25l| of the all-electron linear muffin-tin 
orbital method [26| based on the density functional the- 



ory with local density approximation 27|. The lattice 
constants for Pt and Al used are 3.92 and 4.05 A, respec- 
tively. The basis functions used are s, p, d and / muffin- 
tin orbitals for Pt but s, p, and d muffin-tin orbitals for 
Al [5^. In the self-consistent band structure calculations, 
89 fc-points in the fee irreducible wedge (IW) of the BZ 
were used in the BZ integration. The SHC is evaluated 
by the Kubo formula [1^ . A fine mesh of 60288 fc-points 
on a larger IW (three times the fee IW) is used. These 
correspond to the division of the FX line into 60 seg- 
ments. Comparison with test calculations with 102315 
fc-points (72 divisions of the TX line) for Pt indicate that 
the calculated SHC converges within 1 %. 
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FIG. 1: (color online) (a) Relativistic band structure and (b) 
spin Hall conductivity of fee Pt. The zero energy and the 
dotted line is the Fermi level. The dashed curves in (a) are 
the scalar-relativistic band structure. 

Fig. 1 shows the relativistic band structure of Pt, and 
also the SHC {<Jxy) as a function of Ep- Remarkably, 
the SHC peaks at the true Fermi level (0 eV), with a 
large value of 2200 (fi./e)(f7cm)~^. This gigantic value 
of the SHC is orders of magnitude larger than the cor- 
respondin g va lue in p-type semiconductors Si, Ge, GaAs 
and AlAs [H, Hi]. Furthermore, the calculated SHC in 
simple metal Al is only —17 {H/e){^lcm)~^ , being two 
orders of magnitude smaller than that of Pt. Interest- 
ingly, the SHC in Pt decreases monotonically as the Ep 
is artificially raised and becomes rather small above 3.0 



eV. When the Ep is artificially lowered, the SHC also de- 
creases considerably, and changes its sign at —1.1 eV. As 
the Ep is further lowered, the SHC increases in magni- 
tude again, and becomes peaked at —4.2 eV with a large 
value of —1970 {h/e){Qcm)~^. The SHC decreases again 
when the Ep is further lowered, and finally becomes very 
small below —6.0 eV. Note that the bands below —8.0 eV 
and also above 2.0 eV are predominantly of 5s character 
and the effect of the spin-orbit coupling is negligible. 

We notice that a peak in the SHC appears at the dou- 
ble degeneracies on the L and X points near Ep (0 eV) 
in the scalar-relativistic band structure (i.e., without the 
spin-orbit coupling) while the other peak at —4.2 eV oc- 
curs near the double degeneracies at the L and F points 
(see Fig. [T]). The double degeneracy (bands 5 and 6) at 
L is made mostly (93 %) of dx'z' and d^'z' {z': three- 
fold axis), being consistent with the point group D^d at 
L. The double degeneracy (bands 4 and 5) at X con- 
sists mainly of dx>z> and dy'z' (z': fourfold axis), being 
consistent with the point group D/^h- These double de- 
generacies are lifted by the spin-orbit coupling, with a 
large spin-orbit splittings (^ 0.66, 0.93 eV, respectively). 

One may attribute the large SHC in Pt to these double 
degeneracies. To see this, let us consider the fc-resolved 
contribution to the SHC, i.e.. Berry curvature 17fj(k); 
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where the spin current operator = v}, with spin 

Sz given by Sz = ^P^z {P, ^z- 4x4 Dirac matrices) i28i] . 
/k„ is the Fermi distribution function for the n-th band 
at k. is an analogue of the Berry curvature for the n- 
th band, and it is enhanced when other bands come close 
in energy (i.e. near-degeneracy). Fig. 2(a) shows clearly 
that il^(k) is large only near the L and X points. In- 
terestingly, Berry curvature ri^(k) for the doublet bands 
4 and 5 near the X point are large but have opposite 
signs (Fig. 2(b)). However, because band 5 near the X 
point is unoccupied, only rifi(k) for band 4 contributes 
to the SHC, resulting in the large positive peak in il^(k) 
near the X point (Fig. 2(a)). Fig. 2(c) shows that the 
SHC decreases monotonically as the temperature (T) is 
raised. This rather strong temperature dependence is 
also due to the near-degeneracies since the small energy 
scale is relevant to the SHC there. Nevertheless, the SHC 
a^y = 240(?i/e)(f7cm)"i at T = 300K is still large, and is 
close to the measured value (240) [1]. The SHC for Al at 
4 and 300 K is —17 and — 6(fi/e)(ricm)~^, respectively. 
The former value is similar to the experimental values 
(-27, -34) at 4.2 K [5]. 

In order to study the role of near-degeneracies in more 
detail, we construct two effective Hamiltonians i?(k) for 
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FIG. 2: (color online) (a) Berry curvature r2^(k) at zero 
temperature, and (b) band(n)-decomposed Berry curvature 
57ji(k) along the symmetry lines in the fee Brillouin zone. In 
(b), Q,^(k) for the nth band has been shifted upwards by 
(n — 1) X 500 for clarity. The inset (c) shows the temperature- 
dependence of the spin Hall conductivity a^y. 



1. We flipped the sign of a°'j because the sign of the 
charge of the carriers is opposite from Ref. [sO]. Gfj 
describes the mapping of an area form from the three- 
dimensional k space to the five-dimensional d space. It 
can be regarded as a "solid angle" enclosed by the d 
vector when the wavenumber k runs over the domain 
between the two Fermi surfaces. Hence, it becomes larger 
for smaller d(k) = i(£'u - Ei). The spin operators are 
given by = = T*^/2, and = where 

= ^[r°, r^]. Using these relations, one can calculate 
the SHC cr^ 

points and four L points 

The next issue is whether the contributions from var- 
ious bands cancel or not. From Eq. ([2]), the SHC from 
the X-points and that from the L-points are calculated 
as a function of the Ep, as shown in Fig. [3l Here we 
put a cutoff for the k-integral as 7r/(5a). The integrand 
is dominated by the contribution near the L or the X 
points, and cancellation does not occur when the Fermi 
energy is in the gap. It is analogous to the zero-gap 



from Eq. Q by summing over the three X 



semiconductors rather than GaAs [31|, [Sj] . Thus we can 
identify the peaks at Ep ~ with the peak of the SHC 
in Fig. [U and the enhancement of SHC in Pt is attibuted 
to the near-degeneracies at the L and X points. As is 



the two doubly degenerate bands at X and L points, re- 
spectively. At the X point, by imposing the D^h symme- 
try and the time-reversal symmetry, the effective Hamil- 
tonian with basis |(a;' =F iy')z' 1) and \{x' ± iy')z' j) (z': 
fourfold axis) can be written in terms of 4x4 Clifford F- 
matrices (F^ = t^, — (TzTy, F 
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F^ = a. 
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F^ = T^) as H{k) = e(k) + J2l=i c^a(k)F''. By expanding 
the coefficients da with respect to the wavenumbers k' 
measured from X and L points (k' = k — k;, i — L,X), 
we have constructed the effective Hamiltonian. Fitting 
with the calculated energy bands and wavefunctions, we 
determined the expansion coefficients to fc'^ order. This 
effective model is an even function of k', and is similar 
to the Luttinger model, representing the valence bands 
of cubic semiconductors [BOj , or the valence and conduc- 
tion bands of zero-gap cubic semiconductors [sij near the 
F-point. The previous analysis for the p-type semicon- 
ductors [13] are equally applied. The effective Hamil- 
tonian has the eigenvalues £'i(k) = e(k) — d,(k), and 
Eu(k) — e(k) -I- d(k) for the lower and upper bands, re- 
spectively, where d 



_^ d|, and these bands cor- 



respond to the heavy-hole and light-hole bands, respec- 
tively. From Eq. (35) of Ref. ^3^, the response of a gen- 
eralized spin current (corresponding to F"'') is given by 
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where /^^ and /^i are the Fermi functions of the upper 



and the lower bands, and G' 
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FIG. 3: (color online) Spin Hall conductivity of platinum cal- 
culated from the effective Hamiltonian for (a) the L points 
and (b) the X points, as a function of Ef [33 ]. 



similar to the p-type semiconductors [18|, this intrinsic 
SHE is robust against impurity scattering [3^. To see 
this, we consider dilutely distributed short-ranged impu- 
rities V{y) = V5{r — Vi). It is justified in Pt, because 
screening is prominent compared with semiconductors. 
Then the vertex corrections from the impurity scattering 
for the SHC vanishes in the clean limit from the follow- 
ing reason. Because the effective Hamiltonian satisfies 
iJ(k) — i7(— k), the Green function is an even function 
and the current operator is an odd function of k'. Then 
in calculating the SHC from a correlation function be- 
tween the current jy and the spin current the lad- 
der diagrams from impurities cancel between the inter- 
nal wavenumbers k' and — k' for the current vertex j(k). 
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Thus for short-ranged impurities, the SHC in the clean 
limit is given by the intrinsic value from the bare dia- 
gram without impurity scattering. This justifies our first- 
principle result even in the disordered case. Although it 
may sound trivial, it is not in general; in the Rashba 
model the vertex correction from impurities is relevant, 
and kills the intrinsic SHC even in the clean limit 17 1 . 

We note that i?(k') = H{—'k') results because we re- 
strict ourselves to the even-parity (i.e. d) orbitals. Thus 
even when we include the higher-order terms in k' it holds 
true, and the vertex correction vanishes for short-ranged 
impurities. When k is away from such high-symmetry 
points, the orbitals with odd and even parities are hy- 
bridized, and the SHC will be cancelled to some extent by 
the vertex corrections by impurities. Thus for inversion- 
symmetric systems such as Pt, it is safe to restrict our- 
selves to the high-symmetry points. 

Discussion on the relevance of the present result to the 
experiment on SHE in Pt [3] is in order. At room tem- 
perature the magnitude of a^y ~ 240 {ncm)~^ with the 
conductivity (Txx ^ 10^ ri~^cm~^, corres pon ds to the 
"intrinsic" region in the criterion of Ref. [20[. This is 
consistent with the idea of "resonant" Hall effect since 
the enhanced contribution from the near degeneracies at 
X- and L-points has been confirmed by the present first- 
principles calculation. Hence it is most probable that 
the SHE in Pt at room temperature is due to the intrin- 
sic mechanism calculated in this Letter. On the other 
hand, at the lowest temperature the system enters the 
superclean extrinsic region [S^l, with a^x rising up to 
o'xx ^ 10^ 0~^cm~-'^. Hence at lowest temperature the 
skew scattering becomes very large, and the SHC cannot 
be explained only by the intrinsic mechanism. 

The authors thank National Science Council and 
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CPU time. The work was partly supported by Grant- 
in-Aids (Grant Nos. 15104006, 16076205, 17105002, 
19740177, and 19019004) and NAREGI Nanoscience 
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Note added: After submission of the paper, temper- 
ature dependence of the SHC in Pt was measured to 
be almost constant from T = 300K to T = OK [34 1. 
Though it may look different from our scenario, it is con- 
sistent with it. In 3j|, the conductivity at T = OK 
is cTxx ^ 10^(r2cm)~^, much lower than the above- 
mentioned value. This conductivity corresponds to the 
self-energy of the order of lOmeV, which is comparable 
to room temperature. This implies that the self-energy 
gives a cutoff to the expression of the SHC, and the SHC 
remains constant below room temperature. 
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